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INTRODUCTION 


1.  0 

Subsonic  flow,  pulsed  electric  discharge  lasers  have  propagation  advantages 
and  good  control  over  the  lasing  kinetics.  However,  the  energy  input  to  the 
laser  gas  is  temporally  and  spatially  nonuniform.  This  nonuniform  deposition 
of  energy  into  the  laser  cavity  gas  can  cause  pressure  and  rarefaction  waves 
which,  in  turn,  degrade  the  optical  quality  of  the  cavity.  The  formation  of 
these  waves  thus  places  ultimate  limits  on  the  maximum  pulse  length,  maxi- 
mum repetition  rate,  maximum  output  energy  per  pulse,  and  maximum  output 
power  obtainable  from  pulsed  EDLs. 

The  presence  of  pressure  waves  inside  the  active  volumes  of  electric  discharge 

1 2,  3 

lasers  has  been  previously  recognized  and  attn  _ked.  ’ ’ These  analyses 

used  a Fourier  series  to  represent  the  spatial  distribution  of  the  energy  input 
and  concentrated  on  attenuating  the  fundamental  component.  Thus,  although 
they  were  used  with  good  success  in  designing  the  acoustical  absorbers  for 
repetitively  pulsed  systems,  they  did  not  lend  themselves  well  to  predicting  the 
actual  density  variations  to  be  expected  during  the  lasing  pulse.  Diagnosis  of 
instantaneous  heating  rates  was  rendered  difficult  due  to  the  requirement  of 
solving  for  and  summing  a large  number  of  the  spatial  Fourier  components  of 
the  wave  structure.  In  addition,  they  did  not  attack  the  wave  formed  by  the 
intense  heating  in  the  cathode  fall. 

This  document  summarizes  the  results  of  a combined  analytical/experimental 
effort  at  Northrop  Corporation  under  Contract  N00014-72-C -0043  aimed  at 
providing  improved  and  simplified  means  of  predicting  the  pressure  and 
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rarefaction  wave  structure  found  in  pulsed  electric  discharge  laser  cavities. 


This  analysis  is  addressed  more  to  the  formation  of  waves  during  and  immedi- 


1.  2,  3 

ately  after  the  lasing  pulse  as  opposed  to  the  prior  work  ’ ’ which  addressed 


the  problem  of  attenuation  of  the  waves  after  the  pulse.  In  addition  the 


important  issue  of  heat  transfer  to  the  cathode  from  the  cathode  sheath  has 


been  analytically  attacked  and  its  effect  is  shown. 


For  the  sake  of  generality  of  approach,  a Green's  function  formulation  was 


used  which  lends  itself  well  to  the  approximations  of  uniform  bulk  heating 


and  uniform  cathode  fall  heating.  Simple  expressions  are  thus  derived  which 


can  be  used  to  predict  the  total  density  variations  within  the  EDL  active 


volume.  Using  the  improved  analysis  and  experimentally  determined  density 


profiles,  it  is  feasible  to  predict  instantaneous  heating  rates  during  the  laser 


puise. 
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2.  0 DESCRIPTION  OF  THE  PROBLEM 

Figure  1 shows  a typical  cross  section  of  the  electric  discharge  laser  cavities 
currently  under  study  at  Northrop.  The  optical  axis  is  perpendicular  to  the 
paper  and  the  electric  field  is  parallel  to  it.  An  electron  gan,  at  the  top( 
ionizes  the  laser  gas  with  high  energy  electrons  while  a dc  sustainer  field, 
applied  between  the  anode  and  the  cathode,  provides  the  electrical  power  to 
pump  the  molecules  to  vibrationally  excited  states.  During  the  pumping  and 
optical  extraction  processes,  heating  of  the  gas  occurs.  For  carbon  monoxide 
lasers,  the  dominant  heating  process  is  VV . This  oc  irs  as  a result  of  the 
anharmonic  defect  between  adjacent  vibrational  energy  levels  in  the  gas  and 
the  fact  that  the  vibrational  energy  structure  is  continuously  adjusting  during 
the  pumping  and  optical  extraction  processes. 
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Appearance  of  this  thermal  energy  causes  the  sensible  temperature  of  the  gas 
to  rise  and,  initially,  the  pressure  to  rise.  Since  some  of  the  boundaries  of 
the  cavity  are  either  open  or  porous,  cavity  fluid  leal  s out,  sending  back 
rarefaction  waves  into  the  cavity.  Extremely  high  local  heating  occurs  in  a 
thin  sheath  around  the  cathode.  This  intense  local  heating  causes  an  additional 
wave  structure  to  be  formed  which  sends  a compression  wave  into  the  cavity. 

The  types  of  laser  cavities  studied  in  this  report  have  a porous  anode  through 
which  flow  is  injected  into  the  cavity.  The  sidewalls  are  nonporous  and  guide 
the  fluid  up  towards  the  cathode  screen.  After  the  cathode  screen,  the  flow 
is  directed  sideways  between  the  sidewalls  and  the  electron  gun  to  a large  cross 
sectional  area  manifold  where  it  is  exhausted  by  a pump.  During  typical  pulse 
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time  periods  (t  < t ),  the  exhaust  manifold  may  be  considered  to  be  an  infinite 
volume.  The  same  is  not  true  of  the  channel  in  the  anode  which  feeds  the 
anode's  porous  upper  surface,  however.  In  the  case  of  the  anode,  the  flow 
leaking  from  the  cavity  through  the  anode  screen  will  strike  the  bottom  of  the 
anode  channel  and  reflect  back  towards  the  cavity,  sending  pressure  waves 
into  the  cavity. 

On  the  basis  of  the  above  qualitative  description,  one  would  expect  the  follow- 
ing approximate  distribution  of  waves  in  the  cavity  during  and  after  lasing 
discharge : 

a.  From  the  cathode  area,  a sharp  compression  wave  (due  to 
cathode  fall  heating)  followed  by  a r a refaction  wave  (represent- 
ing leakage  of  cavity  gas  through  the  cathode  structure). 

b.  From  the  anode,  a rarefaction  wave  followed  by  a subsequent 
set  of  compression  waves  caused  by  the  finite  depth  of  the 
channel  underneath  the  anode  screen. 

c.  A final,  sharp  step  compression  wave  on  the  heels  of  the 
cathode  rarefaction  wave s that  would  represent  the  cessation 
of  heating  in  the  cavity  (t  > tp). 

Figure  2 shows  a typical  Mach-Zehnder  inte rfe rog ram  of  a 1 meter  length, 

10  liter  CO  laser  cavity  shortly  after  termination  of  the  lasing  pulse. 
Arj-angement  of  cavity  components  is  the  same  as  in  Figure  1.  As  can  be 
seen,  the  previously  described  wave  patterns  in  a,  , b, , and  c.  , above 
are  clearly  present  in  the  interferogram.  Also  to  be  noted  are  the  relatively 
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Figure  2.  Interferograms  for  electrically 
excited  laser  gas  pumped  at 
(a)  135  W/cm3,  (b)  400  W/cm3, 
and  (c)  1250  W/cm3. 


linear  gradients  and  the  small  percentage  perturbations  in  cavity  density 
(small  number  of  fringe  shifts).  A plot  of  the  cavity  density  gradients  would 
thus  be  expected  to  look  somewhat  like  that  shown  in  Figure  3. 

Because  of  the  relatively  small  percentage  changes  in  cavity  density,  an 
acoustical  treatment  of  the  wave  formation  and  propagation  is  indicated.  The 
linear  density  gradients  imply  relatively  uniform  heating  spatially  and  tem- 
porally. Thus,  the  analytical  formulation  should  be  adapted  to  this  type  of 
energy  input.  These  two  considerations  point  towards  an  acoustical  analysis 
based  on  Green's  functions.  The  Green's  function  approach,  although  quite 
general,  is  well  adapted  to  step  changes  in  heat  rate,  both  spatially  and 
temporally.  The  acoustical  wave  equations,  of  course,  are  linearized. 
Convolution  of  the  Green's  function  solution  with  the  actual  expected  spatial 
heating  profiles  is  thus  valid  within  the  limitations  and  assumptions  inherent 
in  an  acoustical  analysis. 

As  can  be  seen  in  Figure  2,  an  appreciable  and  abrupt  compression  wave 
emanates  from  the  cathode  region.  This  cathode  wave  results  from  the  intense 
heating  in  the  cathode  fall  region  (the  plasma  sheath  which  surrounds  the 
cathode  wires)  and  is  large  enough  to  cause  severe  density  gradients  in  the 
laser  cavity.  Since  it  is  heating  which  causes  the  cathode  wavre  to  be  formed, 
it  is  reasonable  to  assume  that  the  high  temperatures  thus  created  in  the 
cathode  sheath  will  also  transfer  appreciable  amounts  of  thermal  energy  to  the 
cathode  screen  wires.  This  heat  transfer  should  affect  the  magnitude  of  the 
cathode  wave.  Thus,  the  analytical  description  of  the  cathode  wave  structure 
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should  take  into  account,  in  some  manner,  the  effects  of  local  cathode  fall 
heat  transfer. 

The  general  forrr  of  the  analytical  description  of  the  wave  patterns  thus  appears 
to  be  use  of  Green's  functions  and  an  acoustical  analysis.  Heat  transfer  effects 
and  local  heating  in  the  cathode  fall  region  must  be  accounted  for  as  well  as 
porous  boundaries  at  the  top  and  bottom  of  the  laser  cavity.  The  following  two 
sections  outline  the  analysis  and  present  typical  results  for  CO  lasers.  The 
first  ignores  heat  transfer  in  the  cathode  fall  area.  Although  this  leads  to  an 
error  in  the  predicted  results,  the  second  analysis  which  includes  cathode  fall 
heat  transfer  shows  that  this  error  becomes  quite  small  a short  distance  behind 
the  cathode  fall  compression  wave.  The  simpler  zero  heat  transfer  model  can 
thus  be  used  for  interpretation  of  experimental  wave  patterns  including  bulk 
heating  rates.  Because  it  ignores  cathode  heat  transfer,  its  results  are  quite 
easy  to  use. 
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3.0  ANALYSIS  WITHOUT  CATHODE  FALL  HEAT  TRANSFER 
Proceeding  first  with  the  zero  cathode  fall  heat  transfer  case,  we  idealize 
the  actual  cavity  geometry  as  shown  in  Figure  4.  The  cathode  is  assumed 
to  be  a screen  of  very  high  porosity  and  has  the  acoustical  impedance  of  an 
open  channel.  The  anode  is  idealized  as  a porous  plate  of  much  lower 
acoustical  impedance. 


3*  1 Governing  Equations.  To  formulate  the  basic  equations  of  motion,  the 
following  assumptions  are  made: 

a.  Small  amplitude  waves  . 

b.  Negligible  convection  terms. 

c.  Bulk  heating  of  the  plasma  and  local  heating  of  the  cathode  sheath 
are  step  functions  in  time  and  space. 

d.  Perfect  gas  with  constant  specific  heat. 

With  these  assumptions,  one  may  write  the  following  one-dimensional 
equations  of: 

Continuity 


Momentum 


Energy 


BP' 

at 


+ 


0 


0 


(1) 

(2) 

(3) 


and  State 


R[P T ' + TP'] 


(4) 


10 


where  P,  u,  p,  x,  t,  T,  C^,  and  Q have  their  usual  meanings  of  density, 
velocity,  pressure,  distance,  time,  temperature,  specific  heat  at  constant 
volume,  and  heat  input,  respectively.  Primes  refer  to  incremental  changes 
in  quantities,  such  as  velocity  and  pressure,  while  unprimed  quantities 
denote  average  values. 

C ombining  the  above  equations  and  manipulating,  one  obtains  the  final, 
one-dimensional  wave  equation  for  small  pressure  fluctuations: 


( y-  1)  3Q 
2 

c d t 


(5) 


where,  in  addition  to  the  previously  defined  quantities,  7 represents  the 
ratio  of  specific  heat  and  c refers  to  the  speed  of  sound  in  the  undisturbed 
fluid.  An  important  result  of  equation  (5)  is  that  it  is  linear.  Hence,  it  is 
legitimate  to  obtain  separate  solutions  for  bulk  heating  of  the  plasma  and  for 


local  heating  of  the  cathode  sheath  and  combine  the  two  to  obtain  the  complete 
wave  solution  for  the  cavity.  A second  important  result  of  equation  (5)  is  that 
pressure  waves  will  not  be  generated  in  the  cavity  unless  the  second  deriva- 
tive of  heat  addition  with  respect  to  time  is  significant. 


3.2  Boundary  Conditions.  As  a second  order  partial  differential  equation, 
the  wave  equation  above  is  subject  to  two  initial  conditions  in  time  and  two 
boundary  conditions  in  space.  These  conditions,  for  both  the  cathode  heating 
and  for  the  plasma  bulk  heating  solutions,  may  be  written  as  : 


12 


Initial  Conditions  (t  = 0,  all  x): 


= 0 (6) 


3p7/St 

= ( v- 

l)  Q 

(7) 

Boundary  Conditions  (t 

= 0+) 

pu7 

= Kp' 

@ X = 0 

(8) 

p' 

= 0 

@ x = «o 

(9) 

The  two  initial  conditions  refer  to  the  fact  that  the  initial  pressure  fluctuations 
in  the  laser  cavity  are  zero  (6)  and  that  the  initial  rate  of  stored  mechanical 
energy  in  the  gas  is  equal  to  the  initial  rate  of  heat  input  (7).  Customary  con- 
version factors  (J)  between  thermal  and  mechanical  energy  are  omitted. 

The  two  boundary  conditions  refer  to  the  fact  that  the  incremental  flow  through 
the  local  boundary  is  proportional  to  the  incremental  pressure  at  that  boundary 
(8)  and  chat  the  incremental  pressure  changes  far  from  the  flow  boundary  are 
zero  (9).  This  latter  is  equivalent  to  saying  that  the  gas  affected  by  the  rar- 
efaction wave  system  is  of  finite  dimension.  Thus,  for  example,  the  cathode 
heating  is  presumed  to  occur  in  a small  distance,  d,  characteristic  of  the 
cathode  fall  and  this  disturbance  is  not  felt  at  distances,  x,  that  are  far 
removed  from  the  fall  region.  Similarly,  for  the  bulk  plasma  heating  case, 
heating  is  presumed  to  occur  over  a finite  distance,  L,  characteristic  of  the 
cavity  dimension.  At  distances  x > L,  the  analysis  assumes  unheated,  uniform 
uniform  laser  gas  = 0).  For  distances  x»L,  pressure  waves  from  the 
cavity  will  not  have  had  time  to  propagate  or  will  have  attenuated  to  negligible 
value. 


13 


The  factor,  K,  is  a general  acoustic  permeability  of  the  cathode.  If  the  cathode 
has  appreciable  flow  resistance,  then  K will  be  close  to  the  permeability  meas- 


ured under  steady  state  conditions.  As  the  porosity  of  the  cathode  screen 
becomes  very  high,  K will  approach  the  inverse  acoustic  impedance  of  the 
laser  gas  (1  /c). 


Proceeding,  we  may  differentiate  equation  (8)  to  obtain 


du ' 

at 


-K 


at 


which,  when  combined  with  momentum  (2)  yields 

= K (10) 

OX  ot 

as  a second  boundary  condition. 

3,  3 General  Solution.  The  general  solution  for  both  bulk  heating  of  the 
plasma  and  the  local  cathode  fall  heating  consists  of  taking  the  Laplace 
transform  of  wave  equation 

00 

p '(x,  s)  = J c St  p '(x,  t)  dt  (11) 

0 

which  yields  a resulting  one-dimensional  wave  equation  in  distance 


dp'fr,  s)  _ k2  -/ 

dx2 


k p'(x,  s)  = 


O'  - l) 


sQ 


where  k = s/c  and 


6 = J e“St  Q dt 


(12) 


(13) 
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The  solution  is  most  conveniently  accomplished  by  making  a Green's  function 


of  (12)  of  the  form 


2 C'(x/x0)  2 

d - k G(x/Xq)=  6(x  - xQ) 

dx 


where  <5  (x  - x^)  is  the  Dirac  delta  function  and  x^  is  the  position  of  the  unit 
source  in  the  Green's  function.  The  formal  solution  of  (14)  is  divided  into 
two  parts:  x < x^,  and  x > x^;  denoted  G and  G+,  respectively. 


G_  = Cl  d1^  + C2e  ^ (x<  xQ) 


~ + kx 
C1  C 


,A  -kx 
C2C 


(x  >xQ) 


(15b) 


The  Green's  function  is  required  to  satisfy  the  two  continuity  conditions  at 


G-<V 


= G+(xq)  (continuous  through  x^) 


(16a) 


dG  /dx  - dG  /dx  = 1 


(integrated  value  of  the  unit  ( 1 6b) 

im  pulse) 


and,  also,  the  two  boundary  conditions,  (8)  and  (9). 


The  resulting  Green's  functions  are: 


x<  x. 


' [k  (k  + sK)  -][k  cosh(kx)  + sK  sinh(kx)J< 


(17a) 


x > x( 


_ 

G+  = - [k  (k  + sK)  ] [k  CO“h(kx0)  + SK  Sink  "“o']0  ,17!  > 


The  actual  pressure  distribution  is  the  convolution  of  the  Green's  function 
solution  (17)  with  respect  to  the  spatial  distribution  of  heating  (13);  i.e., 


(x,  s)  = - |^(  ^ 2~— j j s 6 G(x/xq)  dxf 


Inversion  of  (18)  yields  the  final  desired  solution  of  incremental  cavity  pressure 
as  a function  of  time  and  distance. 


3.  1 Cathode  Wave  Solution.  Proceeding  first  with  the  cathode  wave,  we 


represent  the  spatial  heating  profile  by  the  expressions 


0(x,  t) 


constant 


x < d 


(19a) 


Q(x,  t) 


x > d 


'19b) 


and  the  temporal  heating  profile  by  the  expressions 


constant, 


t < t 


(20a) 


t > t 


(20b) 


where  t refers  to  the  time  that  electrical  discharge  actually  occurs  in  the 
P 

laser  cavity. 


Substituting  equations  (19)  and  (20)  into  (18),  integrating  and  inverting,  yields 
the  final  expression  for  the  pressure  wave  system  generated  by  the  cathode 


S ■ i »(>-“')[•  - * -Sfe  "(••“t  ■ f1] 


-T7ET  “KIM] 


where,  in  addition  to  the  previously  defined  quantities,  U(  ) refers  to  the 


unit  step  function. 
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The  above  equation  may  be  simplified  in  notation  by  denoting 


* 

t = 

t/lP 

(22a) 

4 

X = 

x/(t  c) 
P 

(22b) 

¥ 

P 

P 'I? 

(22c) 

• 4 

Q 

c 

(Q  t )/p 

C p 

(22d) 

* 

d 

d/(tpC) 

(22e ) 

with  the  aid  of  which,  equation  21  may  be  rewritten  in  the  following  nondimen- 
sional  form : 


4 4 4 


p*(x*,t*)  = (V-DQcjfMdV)  U(t  +d  -X  ) 


- (t*-d  -x  ) U(t  -d  -x  ) - (t‘  -x“  ) U(t^-x  ))  (23) 

2(1+Kc)  V 1+Kc  J 


An  interesting  and  important  result  may  be  obtained  by  letting  the  thickness, 
d,  of  the  cathode  fall  region  approach  zero.  Equation  (23)  is  then  reduced  to 

the  form  . * * 

Qc  d 

Lim  p = (y-  1\TTKc)  (23 

d-0  1 

In  this  limiting  case,  the  amplitude  of  the  cathode  wave  depends  only  upon 
the  total  power  put  into  the  cathode  fall  region  and  not  upon  the  details  of 
how  deep  the  cathode  sheath  is  or  upon  the  local  heating  rate.  Since  the 
cathode  sheath  will  be  very  small  compared  to  the  dimensions  of  the  laser 
cavity  and  the  overall  wave  structure  inside  it,  this  limiting  case,  with  its 
attendant  simplified  analytical  results,  becomes  a reasonable  description  of 

the  cathode  wave  structure. 


.1 


1 


■J8 
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3.  5 Bulk  Heating  of  the  Plasma.  The  relatively  straight  and  uniformly 
spaced  fringes  of  the  central  areas  of  the  interferograms  in  Figure  2 suggest 
uniform  heating  of  the  plasma  in  the  laser's  active  volume.  Thus,  for  the 
sake  of  ease  of  analysis,  we  have  assumed  that  the  bulk  heating  of  the  plasma 
may  be  represented  by: 

Q = 0 for  t < 0 (24a) 

Q = Qh  = constant  for  t > 0 (24b) 

Equation  (24)  may  be  substituted  directly  into  (18)  to  obtain 

(V-D  Qx 


(x,  s)=  LimJ  - * s Q G(x/xQ)  dxQ=  2 [T  1+Kc 

d— ®°  o c s 


-sx 
c 


|l-—  1 

L l+Kc  J 


(25) 


the  inverse  of  which  yields  the  following  expression  for  the  bulk  heating 
induced  expansion  wave  which  emanates  from  the  cathode 

Pc  <**•  11  * '"-"“h  ['  ■ (t  ••  x ' U(t" x ’] 


(26) 


An  expansion  wave  also  emanates  from  the  anode  due  to  the  anode's  finite 
porosity.  Denoting  the  effective  acoustic  permeability  of  the  anode  by  Ka  and 


substituting  L-x  in  place  of  x above,  we  may  write  the  nondimensional  anode 


wave  as 


* * * • * r * 

pa(x  , t ) = (?-  1)QU  [t  - f 


K c 
a 


+ K c 
a 


(t  "+  x - L )U(t  + x - L )]  (27) 


Since  the  differential  equations  are  linear,  the  complete  pressure  wave 
system  in  the  laser  cavity  may  be  obtained  by  adding  the  individual  solutions 
for  the  cathode  fall  wave,  pf'i  the  bulk  heating  induced  cathode  wave,  p , and 

I C 


i 


i 
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the  bulk  heating  induced  anode  wave,  p . Combining  equations  (23),  (26), 
and  (27)  yields  the  final,  complete  solution  for  the  pressure  fluctuations 
inside  a pulsed  electric  discharge  laser's  cavity  during  and  immediately 
after  discharge: 


K c 
a 


(X  , t ) = (y-  l)Qh  (t  - (t  -X  ) U(t"-x " )-  (t+x-L)  U(t  +x  -L  ) 

1 * a 


-jj(t  +d  -x  ) U(t  -Ld  -x  ) - “r777Tv(t  -d'-x  ) U(t  -d  -x  ) 


Q 

n 

Kc 

1+Kc 


2(  1 f Kc) 


(t  -x  ) U(t 


■-•»]) 


(28) 


3.6  Density  Fluctuations  in  the  Laser  Cavity.  It  is  the  density  fluctuations 
in  the  cavity  which  are  of  primary  concern  since  it  is  these,  not  pressure 
per  se,  which  affect  the  laser  beam.  Fortunately,  within  the  limits  of  small 
perturbation  theory,  the  density  fluctuations  may  be  simply  related  to  the 
previously  derived  pressure  fluctuations.  Using  the  first  law  for  heat  addition 
into  the  gas: 

cdT  = +4^  (29) 

V P p c. 

and  the  equation  of  state,  the  density  waves  will  be  of  the  form 

dp  = lE  . (ZlJJ  dQ  (30) 

2 2 h 

c c 


o r, 


using  primed  quantities  instead  of  differentials 


I 


Performing  the  above  substitution,  equation  (28)  may  be  rewritten  in  terms 
of  the  ordered  density  profiles  to  be  expected  in  the  laser  cavity: 


[p  -«v-  ■»<  1 ] 

Kc 


2^1 6 
y n 


K c 


l + Kc 


(t  -x  ) U(t  -x  ) - 


1+K  c 
a 


(t  +x  -L  )U(t  +x  -I.  ) 


+ — [|(t  +d  -x  ) U(t  +d  -x"  ) - 2(1+^) (t  _d  'X  ) U(t  "d  'X  } 
Q. 


(1-Kc) 


- TTKc(t  -x  )U(t 


-*>]) 


(31) 


Although  somewhat  long,  equation  (31)  is  analytic  and  easily  evaluated.  A 
typical  set  of  results  is  shown  in  Figure  5 in  which  normalized  density,  P , 
is  plotted  against  normalized  distance  x , for  various  normalized  times, 
t . The  cathode  is  assumed  to  be  completely  porous.  Hence,  K = 1/c. 
Other  assumed  parameters  are  listed  on  the  figure. 


II 


The  cathode  wave  structure  propagates  from  the  left.  It  is  composed  of  both 
the  cat  Kde  sheath  heating  and  the  rarefaction  wave  which  results  from  plasma 
bulk  heating.  The  step  increase  in  density  results  from  the  assumption  that 
heating  in  the  cathode  sheath  is  uniformly  distributed.  Its  amplitude  meas- 
ures the  local  energy  flux,  O^d,  into  the  cathode  sheath.  The  linear  slope 
following  this  step  rise  indicates  that  bulk  heating  of  the  plasma  is  uniform 
and  that  the  acoustic  impedance  of  the  cathode  is  constant  (that  of  the  laser 
gas  itself).  At  time  periods  exceeding  the  plasma  discharge  (t/t^  > 1.0),  a 
step  decrease  in  the  density  appears.  This  step  decrease  is  a negative  wave 
produced  as  the  local  plasma  sheath  around  the  cathode  disappears. 
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Cathode  Waves 

Vs 

■ 14  kV 

!s 

= 1000  Amp 
■ 370W/cm2 

°cd 

T0 

- 300°K 

Po 

c„ 

a 400Torr 
• 3.  5 x 104  cm 

Figure  5.  Theoretical  density  fluctuations  in  an  electric  discharge 
laser  cavity.  No  heat  transfer  to  cathode. 


To  be  noted  are  the  different  slopes  of  the  cathode  and  anode  wave  structures. 
These  different  density  slopes  reflect  the  difference  between  the  assumed 
permeabilities  for  the  cathode  and  anode  (1:  1/23).  As  the  bulk  heating  rate 
is  increased,  these  slopes  will  also  change  as  shown  in  Figure  6.  Their 
ratio,  however,  will  remain  the  same.  The  initial  height  of  the  cathode  wave 
will,  of  course,  remain  constant  as  long  as  the  assumed  power  input  to  the 
cathode  fall  remains  constant. 

Figure  7 shows  a comparison  between  actual  density  profiles  as  measured  in 
a Northrop  10  liter  EDL  and  those  predicted  by  this  zero  heat  transfer  analysis 
(A  heating  rate  of  15%  of  imput  power  was  assumed.  ) The  slopes  of  the  wave- 
fronts  and  their  propagation  into  the  cavity  compare,  in  the  case  shown,  within 
a few  percent.  The  main  region  of  uncertainty  appears  to  be  in  the  front  of 
the  cathode  wave  which  the  interferogram  is  not  able  to  resolve. 

This  agreement  between  experiment  and  theory  is  rewarding  since,  as  previ- 
ously stated,  large  amounts  of  heat  transfer  to  the  cathode  screen  should  occur 
which  the  present  analysis  neglects.  Thus,  this  simplified  analysis  should  be 
valid  for  interpretation  of  the  inte rfe rograms  of  EDL  discharges  and  resultant 
predictions  of  instantaneous  bulk  heating  rates.  It  is,  however,  important  to 
know  the  errors  and  approximations  introduced  by  neglecting  heat  transfer  to 
the  cathode  screen  wire.  The  following  section  briefly  outlines  the  complete 
analysis  including  cathode  sere  *n  heat  transfer  and  points  out  the  essential 
differences  which  result. 


Figure  6.  Illustration  of  the  effect  of  changing  the  bulk  heating 
ratio.  Zero  heat  transfer  to  the  cathode. 


DENSITY  VARIATIONS  PRODUCED  BY 
CATHODE  AND  ANODE  WAVES 


Figure  7.  Comparison  of  the  experimental  and  theoretically  predicted  density  variations 
in  a CO  electric  discharge  laser.  Zero  heat  transfer  to  the  cathode  assumed. 


4.  0 ANALYSIS  WITH  CATHODE  FALL  HEAT  TRANSFER 

As  stated  earlier,  high  local  heating  rates  and  possible  high  temperatures 

are  to  be  expected  in  the  plasma  sheath  which  surrounds  the  cathode  wires. 

As  an  example,  one  may  assume  a typical  cathode  fall  of  roughly 

V = 370V 

c 

a sheath  thickness  of 


d = 0.  1 mm  = 0.  01  cm 


and  an  average  current  density  of 

I/A  = t.OA/cm2 

Actual  local  current  densities  on  the  individual  cathode  screen  wires  will, 
of  course,  be  much  higher  since  their  area  is  not  spread  uniformly  over  the 
entire  discharge.  Assuming  nitrogen/CO  for  the  gas  mixture,  a 100  micro- 
second pulse,  and  an  average  pressure  of  400  Torr,  one  then  calculates  a 
typical,  average  sheath  heating  rate 


6 = (Q  d)/d  = 370/(0.  01)  = 37,  OOOW/cm3 

c c 

and  a theoretical  maximum  temperature  rise,  assuming  a 100  microsecond 
pulse,  of 


T 

max 


Q /(PC  ) = 

C V 


37, 000 

(0.  6 x 10'3)(0.  718) 


= 8,  550°C 


where  P is  the  nominal  gas  density  and  is  the  gas  specific  heat  at  constant 
volume. 


This  high  temperature  rise  cannot  exist  in  practice.  Using  nominal  gas 
thermal  conductivities,  half  the  cathode  fall  thickness  as  a typical  distance, 
one  calculates,  at  the  end  of  a 100  microsecond  discharge,  heat  transfer 
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rates  to  the  cathode  that  are  in  excess  of  the  input  power  to  the  cathode  fall 
itself.  It  is  thus  apparent  that  heat  transfer  to  the  cathode  wires  is  signifi- 


cant and  will  substantially  alter  the  details  of  the  compression  wave  generated 


by  the  cathode  fall. 


Qualitatively,  one  may  a priori  estimate  the  influence  of  cathode  sheatn  heat 


transfer.  At  the  beginning  of  the  discharge,  when  temperatures  are  low,  heat 


transfer  effects  should  be  insignificant.  Thus,  the  initial  height  of  the  com- 


pression wave  emanating  from  the  cathode  should  remain  unchanged.  At  the 


end  of  the  laser  discharge,  when  adiabatic  temperatures  would  be  very  high, 


heat  transfer  effects  should  dominate  and  the  cathode  fall  compression  wave 


should  be  considerably  weakened.  Thus,  one  would  expect  the  composite 


density  profile,  including  the  effects  of  heat  transfer,  to  look  much  like  that 


sho-vn  in  Figures  5 and  6 except  that  the  slope  behind  the  advancing  compres 


sion  wave  on  the  left  hand  side  of  the  figure  would  now  not  be  constant. 


Hopefully,  the  change  in  slope  should  die  out  relatively  quickly  and  one  would 


be  left  with  a slope  characteristic  of  the  impedance  of  the  gas  on  the  E-gun 


side  of  the  cathode.  Results  at  that  point  would  be  identical  with  those  shown 


in  Figures  5 and  6. 


To  analyze  the  cathode  fall  produced  compression  wave  with  heat  transfer, 


we  must  write  the  full,  linearized  energy  equation 


dT ' p c>u ' 6 ^ d T \ 

p — + tv  — ■'  C~v  + cv  l.x2  > 
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This,  when  combined  with  the  equations  of  continuity,  momentum,  and  state, 
(equations  (1),  (2),  and  (4),  respectively),  yields  a modified  pressure  wave 
equation  which  accounts  for  heat  transfer: 

2 . . 2 

& x 


4-'-  o--i)|f+(r-i) 


dt 


o x 


JL  (±JL!  \ 

St  L 2 ) 


(33) 


As  can  be  seen,  equation  (33)  is  identical  to  the  adiabatic  wave  equation 
except  for  a temperature  gradient  term  on  the  right  hand  side.  Unfortunately, 
this  temperature  term  is  not  decoupled  from  the  pressure  terms  on  the  left 
hand  side  of  the  equation.  The  simplified  Green's  function  approach  is  no 
longer  appropriate. 


This  dilemma  may  be  resolved  by  assuming  that  heat  transfer  effects  are 
important  only  in  the  cathode  sheath  and  that  this  distance  is  very  small  com- 
pared to  the  characteristic  dimensions  of  the  laser  cavity  and  the  resultant 
gross  wave  structure.  By  so  doing,  one  is  then  able  to  obtain  a separate 
approximate  solution  for  the  conduction  dominated  heat  transfer  temperature 
profiles  in  the  cathode  sheath  and  substitute  these  profiles  into  the  tempera- 
ture term  on  the  right  hand  side  of  (33).  Since  the  obtained  temperature 
profiles  will  be  a function  of  the  input  heating  rate  only,  the  wave  equation 
again  becomes  explicit  in  pressure  and  may  be  solved  in  the  same  manner 
as  used  for  the  adiabatic  case. 

Proceeding  in  this  manner,  we  analyze  the  temperatures  within  the  cathode 
sheath  by  making  the  following  assumptions: 
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a.  The  surface  temperature  of  the  cathode  wires  does  not  change 
appreciably  during  the  discharge  time  period;  and 

b.  Convection  terms  can  be  neglected  in  the  sheath  region. 

The  first  assumption  may  be  justified  by  making  a simple,  one-dimensional 
slab  heating  calculation.  Because  of  its  much  higher  density,  the  tempera- 
ture rise  on  the  cathode  wires'  surface  required  to  absorb  the  cathode  fall 
heat  transfer  is  quite  minimal  - in  the  order  of  a few  degrees  centigrade. 
The  second  assumption  is  made  primarily  to  ease  the  analysis.  The  resul- 
tant analytical  solution  is,  however,  consistent  with  this  assumption  and 
indicates  that  heat  transfer  effects  outweigh  convective  work  in  the  cathode 
fall  plasma  sheath. 

With  the  aid  of  the  above  assumptions,  the  energy  equation  for  the  cathode 


fall  may  be  written 


PC 


dT_' 
v d t 


. aZT' 

X ~l  + 

d x 


Q 


(34) 


This  energy  equation  differs  from  that  previously  used  (equation  (3))  by 
neglecting  convective  mechanical  work  (pfdu/dxJ/C^)  and  substituting,  in 
its  place,  conductive  heat  transport.  By  so  doing,  the  energy  equation 
becomes  decoupled  from  pressure  terms  and  hence  may  be  solved  separately 
from  the  wave  equation. 


The  solution  of  (34)  is  most  conveniently  accomplished  with  the  aid  of  Green's 
functions.  Boundary  conditions  for  the  solution  are: 
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(?  x = 0,  all  t 


@ t = 0,  all  x 


@ x — oe 


Taking  the  Laplace  transform  of  (34)  we  obtain 

a2?  „ 2 t S_ 

— 2 - ^ T - - x 

d x 

> 

where  A is  the  gas  thermal  conductivity  and  ft  is  the  expression 

PCs 

~2  v 


The  resulting  Green's  function  is  then  found  to  be 


= - Ja/s  exp  [-XqA/o/s]  s inh( yS/Q  x)  x < Xq 

Ci*  = -yjajs  exp[-x  / /a/s]  sinh(  ^s/a  xQ)  x > xQ 

and  the  resulting  temperature  profile  is  given  by  the  expressions 


x < d 


x > d 


T (x,  s)  = - a Jexp[-  y/s/a  x]  - 1+-|  exp[ys/a  (x-d)] 
- J exp[->/s/a  (x+d)j| 

Qc  r 

T+  (x,  s ) = - — y Q exp[-  y/s/a  x]  jcosh  (y/s/a  d)  - lj 


(38a) 


(38b) 


As  can  be  seen  above,  the  conduction  dominated  temperature  profiles  in  the 
cathode  sheath  are  a function,  3 V -ly,  of  the  cathode  fall  heating  rate,  Qf  . 
Furthermore,  they  die  out  very  rapidly  - in  a distance  characteristic  of  the 
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cathode  fall  depth,  d.  Thus,  using  the  previously  described  artifice  of  sub- 
stituting the  temperature  profiles  obtained  above  into  the  wave  equation  (33), 
one  obtains  an  explicit  wave  equation  with  the  same  Green's  functions  as 
obtained  for  the  adiabatic  solution.  The  resultant  approximate  solution  for 
the  cathode  fall  compression  wave  with  heat  transfer  effects  then  can  be 
obtained  by  convolving  the  adiabatic  Green's  functions  with  the  forcing  term 
on  the  right  hand  side  of  (33): 


?(x,  s)  = - 


[Ll 


nf(-  5- 


2 

c 0 


Xs  — y|g  (x/xq)  dx 
d x 


0 


(39) 


Vroceeding  with  the  indicated  integration,  we  obtain  for  distances  x > d, 

2 

, f d T , 

G — dxQ  + As J G — dx 

d v 

0 


d d 9 ^ x 2 ■ 


R G+d-V  KsIC 


dx 


0 


« h2t+ 

+ Xs  / G‘  dx 

v dx^ 


(40) 


For  simplicity  we  assume  that 


as 


« l 


(41) 


with  the  aid  of  which  equation  (40)  may  be  integrated  to  obtain 

o 1 e_kxri+Kc  kd  (1-Kc)  -kd  T,  1 

— c — — — e - e - Kc 

(V-l)Q  (1+Kc)  2 12  2 J 

c s 


uTko  [(^T7r)cxp[-sx/c  • ^ad]  - ^372cxpl-sx/c)J 

~~2  ^exp[  ~JS  (x-d)]+  exp[  -y/s  /a  (x+d)]  - 2 exp[  - >/s  j ax  ] 


2s  c 


- exP  [ 


-s(x-d)/c]  +(Tf^)oxpl-s<x+d)/c'+(7)fl^r)cxP<  Tl]  (42) 


30 


The  inverse  Laplace  transform  of  (42)  may  be  separated  into  two  parts:  the 
previously  obtained  adiabatic  wave  solution  (equation  (28))  and  a solution 
which  includes  the  effect  of  cathode  fall  heat  transfer.  This  latter  solution 
is  denoted  p^  and  is  given  below  in  equation  (43). 


where 

, * Q , X 

° „ — and  a = — (■ 

c tp  V 

and  U(  ) is  the  previously  described  unit  step  function. 

To  illustrate  the  effects  of  cathode  sheath  heat  transfer,  the  complete 
pressure  wave  equation [(43  ) plus  the  adiabatic  solution  ( 28 )] is  plotted  in 
Figure  8.  Figure  9 plots  the  corresponding  density  profile.  Several 
important  results  stand  out.  First,  looking  at  the  zero  cathode  porosity 
case  (K  = 0),  it  is  seen  that  the  initial  height  of  the  cathode  wave  is 
unattenuated  by  heat  transfer.  This  is  an  expected  result  since,  at  t = 0, 
no  gas  heating  has  occurred  and,  hence,  no  heat  transfer  could  take  place. 


-A-  = 0. 587  cmwsec 
PCV 

t/tp  - 0. 4 

All  other  parameters  are 
the  same  as  Figure  5. 


K a 0 


without  Heat  Transfer 

i 
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J''' 


Figure  8.  Typical  cathode  fall  produced  pressure  profiles 
with  and  vvithout  heat  transfer  to  the  cathode. 
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Figure  9.  Typical  cathode  fall  produced  density  profiles 
with  and  without  the  effects  of  heat  transfer. 


As  time  progresses  (x/L  becomes  less  in  the  figure),  the  magnitude  of  the 
wave  decreases  since  more  of  the  fall's  energy  is  going  into  heating  the 
cathode  screen  wires  as  opposed  to  mechanical  work  on  the  surrounding 
laser  gas.  After  some  time  has  passed  (x/L  = 0.05  in  the  figure),  a steady 
state  has  been  established.  The  rates  of  cathode  fall  power  input  into 
cathode  wire  heating  and  into  mechanical  work  on  the  gas  are  constant 
percentages  and  remain  so  for  the  duration  of  the  discharge  time  period. 

For  the  high  porosity  cathode  case  (K  = 1/c),  the  same  qualitative  effects 
occur  except  that  the  open  area  behind  the  cathode  permits  gas  to  Tow  from 
behind  the  compression  wave,  thus  relieving  the  pressuie  and  weakening 

the  wave. 

The  important  point  to  be  drawn  from  the  two  figures  (8  and  9)  is  that  the 
slopes  of  the  pressure  behind  the  cathode  wave  are  essentially  identical  for 
both  the  adiabatic  and  the  finite  cathode  heat  transfer  cases  after  a short 
distance  behind  the  compression  wave  front.  In  the  case  snown,  this  distance 
is  roughly  5%  of  the  total  cavity  length.  The  percentage  of  the  total  wave 
-.tructure  affected  by  cathode  fall  heat  transfer  is  thus  quite  small  and  the 
majority  of  the  wave  pattern  pressure  and  density  gradients  are  identical 
for  both  the  adiabatic  and  the  finite  heat  transfer  cases.  This  fortunate 
concurrence  enables  one  to  interpret  the  interferograms  made  of  a pulsed 
EDL's  cavity  with  the  aid  of  the  simpler  adiabatic  density  wave  analysis. 
Typical  errors  using  the  adiabatic  analysis  should  be  no  greater  than  1%  or 
2%.  Inaccuracies  incurred  as  a result  of  using  this  approach  are  thus  much 
smaller  than  those  occurring  experimentally. 
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DETERMINATION  OF  BULK  HEATING  OF  THE  LASER  PLASMA 


5.  0 


Using  the  previously  developed  predictions  of  the  density  profiles  in  the 


laser  cavity,  prediction  of  the  bulk  heating  of  the  laser  plasma  from  experi- 
mental interferograms  becomes  a relatively  straightforward  task.  The  formal 
approach  would  be  to  differentiate  every  term  of  the  adiabatic  wave  equation 
(equation  (31))  with  respect  to  x.  Use  of  the  adiabatic  wave  equation  is 


justified  by  the  previous  section  in  which  it  was  shown  that  cathode  sheath 
heat  transfer  effects  on  the  wave  structure  die  out  a short  distance  behind 


the  compression  wavefront. 


I 


. 

[ 

k 


Fortunately,  it  is  not  necessary  to  differentiate  every  term  of  (31).  The 

equation  is  split  into  two  parts:  The  first  two  terms  describe  the  wave 

patterns  caused  by  bulk  heating.  The  second  set  of  terms,  preceded  by  the 

ratio,  (Q  /Q.  ),  describes  the  wave  pattern  produced  by  the  cath  .1 

c h 

heating.  This  second  set  of  terms,  if  differentiated,  cancel  ide  .cally. 
This  amounts  to  saying  that  the  cathode  fall  heating  causes  a steep  change 
in  density  profile  and  does  not  contribute  to  density  gradients  behind  the 
compression  wave  front. 


Performing  the  indicated  differentiation,  one  obtains,  for  the  cathode 
rarefaction  wave: 


dp  / dx ) 


cathode 


O'-  1)  Qh  (Kc) 
y cP  (i  + Kc) 


(45) 


! 


i 
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and,  for  the  anode  rarefaction  wave: 


dp  /dx) 


( y - i) q.  <k  c) 

h a 


anode 


cp  (1  + K c) 
a 


(46) 


Bulk  heating  of  the  plasma  is  usually  defined  in  terms  of  a percentage  heat- 
ing rate,  *7  : 


rj  = 0,/E 

h in 


(47 


Substituting  this  into  (45)  and  (46)  and  noting  that  the  cathode  screen  perme- 
ability is  essentially  that  of  the  laser  gas,  1/c,  we  may  write 

yep  (1  + Kac)  (dP  /dx)anode  zy cp  (dP  /dx)cathode 


T)  = 


(Y-  1)  k c E. 

a in 


(y-  1)  E. 


in 


The  effective  heating  rate,  H , may  thus  be  checked  by  both  the  cathode  and 
the  anode  wave  patterns  providing  that  the  acoustic  permeability  of  the  anode 
is  sufficiently  well  known. 


The  cavity  density  gradient  is  most  conveniently  obtained  from  a Mach- 
Zehnder  interferogram.  With  the  aid  of  the  interferogram,  the  density 
gradient  is  simply  expressed  as: 

( P / p ) A dN/dx 

dP/dx  = 5 fz W 

where  P is  the  density  of  the  gas  at  standard  conditions,  p is  the  gas 
s o 

density  at  test  conditions,  L is  the  optical  path  length,  Q is  the  Gladstone- 
Dale  constant,  A is  the  wavelength  of  the  light  used,  and  dN/dx  is  the  fringe 


gradient  on  the  interferogram. 


A Mach-Zehnder  interferometer  was  set  up  to  measure  density  waves  in 
the  Northrop  10  liter  CO  electric  discharge  laser.  The  setup  is  schemati- 
cally shown  in  Figure  10.  A ruby  laser  was  used  for  making  single  and 
double  pulsed,  time  resolved  interferogram s and  a helium-neon  laser, 
mouiiced  on  the  ruby  laser's  optical  axis,  was  used  for  system  alignment. 
Firing  of  the  ruby  laser  was  correlated  with  firing  of  the  CO  laser  by  means 
of  a fast  rise  time  transducer  fed  with  a fraction  of  the  ruby  laser's  output. 


This  trace  was  displayed  in  parallel  with  the  sustainer  current'  of  the  CO 

laser  on  a dual  trace,  high  speed  oscilloscope.  Details  of  the  experiment 

are  described  in  Reference  6.  Details  of  the  Northrop  10  liter  laser  are 

5 

reported  in  a previous  semiannual  technical  report. 

Discharge  measurements  were  made  on  pure  N^  CO/Ar  mixes,  and  CO/N2 
mixes.  No  optical  extraction  was  made.  Figure  11  plots  observi  d fringe 
gradient  behind  the  cathode  wave  as  a function  of  input  power  for  various 
combinations  of  sustainer  voltages  and  gas  mixtures.  Figure  12  plots  the 
derived  heating  rate,  T], 

From  Figure  12,  several  important  conclusions  can  be  drawn  regarding 
bulk  heating  of  the  laser  discharge,  at  least  in  the  absence  of  optical 
extraction. 

1.  At  high  sustainer  voltages,  the  bulk  heating  rate  of  both  N^ 


and  the  N^/CO  mixture  depended  very  little  upon  input  power 
for  pumping  rates  above  100  W/cmJ. 
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Experimentally  determined  fringe  gradients  from  interferograms  of  the 
Northrop  10  liter  CO  electric  discharge  laser.  T:>  = 100  Torr,  T = 80°K 


2.  At  low  sustainer  voltages,  heating  rate  does  increase  for 


increased  pumping  rate. 

3.  Over  the  limited  range  of  similar  voltages  tested,  the  heating 
rate  for  N^/CO  does  not  appear  to  be  appreciably  different 
than  that  for  alone. 

4.  Similarly,  for  the  limited  range  and  low  sustainer  voltage 
tested,  heating  rates  for  CO/Ar  and  CO/N2  appear  to  be 
approximately  the  same. 

5.  Heating  rate  increases  as  sustainer  voltage  is  lowered  for 
the  CO/N^  mixes. 

6.  In  general,  the  heating  rates  for  CO  lasing  gas  mixtures 
appear  to  be  quite  low.  In  all  cases,  they  are  below  20% 
and,  for  maximum  energy  sustainer  voltages  (12-13  kV), 
heating  rates  are  below  10%. 


This  last  conclusion  is  particularly  rewarding  since  it  provides  direct  con- 
firmation of  the  low  heating  rates  predicted  by  the  Northrop  CO  laser 
Ivinetics  code  and  thus  confirms  the  potential  high  efficiency  of  CO  lasers 
in  general. 
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